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CHAPTER 1 


INTRODUCTION 

The double oxidation techriiaue usina 0-18 as a 
tracer has been successfully used to identify the 
transport mechanisms in growing oxide scales. The 
detailed mechanism of oxidation is of more than 
academic interest as it determines growth stresses 
and oxide spallation. The aim of this study was to 
develop analysis technioues to obtain Quantitative 
diffusion data from the tracer concentration 
profiles in growing oxide scales. An oxide can grow 
inward at the metal-oxide interface by inward 
oxygen diffusion* it can grow outward at the gas 
oxide interface by outward cation diffusion or it 
can grow by a combination of both the mechanisms. 
The problem of determining the fraction of inward 
and outward growth of the scale from the tracer 
concentration profile was solved. The process of 
tracer diffusion in a growing polycrystalline scale 
was theoretically modelled. 
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The analysis technioues developed were then 
applied to a FeCrAl-Zr alloy by using SIMS to 
obtain the tracer profile* 
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CHAPTER 2 


LITERATURE SURVEY 

2.1 Mathematics of diffusion 

2.1.1 Diffusion 

Atoms in a solid are in constant motion about 
their lattice sites. Occasionally* a few atoms 
sain sufficient energy to Jump into neighbouring 
sites. This random movement of atoms is called 
diffusion. Self diffusion is the phenomenon of 
migration of atoms in the absence of concentration 
gradients. Self diffusion can be studied by 
introducing an isotopic tracer into the system 
such that the diffusing atoms can be identified. 
Strictly speaking the diffusion coefficients of 
isotopes are related by the eauationJ 



Where M and M are the masses and D and D are the 
12 12 
diffusion coefficients of the isotopes. However* 

generally the isotopic masses are so close to 

each other that the diffusion coefficients are 

assumed to be the same. When 0-18 is used as a 

tracer to study oxygen diffusion* the difference 
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is about 6V. 


2.1.2 Fick's laws 


Matter 

flows in a 

manner to 

decrease 

its 

concentration gradient 

or more 

precisely 

its 

chemical 

potential 

gradient. 

The 

f lux 


r> 


density. J (ks/m sec) » of the species i in 3 
i 

dimensions is given by! 

J i " _D i P VC i ...(2.2) 

o 

Where D is the diffusion coefficient <m /sec). 

3 

P is the density (ks/m ) and C i is the mass 
fractional concentration gradient <m _1 ). This is 
Pick's first law. In 1-dimension Equation 2.2 
becomes ! 


„ 3C. 

-D. p i 

1 3x 


..(2.3) 


From Eauation 2.3 and using the principle of 
mass conservation C13 it can be shown that! 

8 C i 9 (D, 3C i \ 

- \ i 1 ...(2.4) 

3 t 3x 3x 

where t is the time* This is Fick's second law. 


If D is independent of position t Eauation 2.4 
i 

becomes ! 


3C . 

-5 r ■ D i 


3c i 

"^ 2 


.(2.5) 


2.1.3 Diffusion mechanisms 
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Diffusion generally occurs by a vacancy or an 
interstitial mechanism. An atom adjacent to a 
vacancy can exchange sites with the vacancy and 
thus migrate through the crystal. Some atoms 
occupy positions in the interstices of the 
crystal lattice. They can move either by Jumping 
from an interstitial site to another (interstitial 
mechanism) or by displacing a lattice atom into an 
interstitial position and taking up the lattice 
site itself ( intersti tialcy mechanism). 

The predominant point defects in an ionic 
crystal are C21J 

i) Schottky disordei — Stoichiometric proportion of 
cations and anion vacancies. 

ii) Frenkel disorder- Eaual number of cation 
vacancies and interstitials. 

iii) Anti Frenkel disorder- Eoual number of anion 
interstitials and vacancies. 

iv) Anti Schottky disorder- Stoichiometric 
proportion of anion and cation interstitials. 
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It should be noted that the experimental 
diffusion coefficient can be represented by an 
Arrhenius type eciuation* 

D » D o exp (- ...(2.6) 

r> 

where D is the preexponential term (m /sec) and Q 
o 

is the activation energy (kJ/mole). 

2.2 Grain boundary diffusion 

2*2* 1* Introduction 

Enhanced grain boundary diffusion has been 
observed in several oxides C3-4H. The observed 
tracer concentration profile in a polycrystsl 1 ine 
sample is due to both fast diffusion along the 
grain boundaries and slower diffuson into the 
grains from the surface and the grain boundaries. 
Some of the models Proposed to mathematically 
treat the combined grain boundary and volume 
diffusion problem are discussed in the next 
section. 

2*2.2 Mathematical treatment 
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Fisher E53 modelled a grain boundary as a 


thin slab of thickness w extending into a semi 

infinite crystal (Figure 1> The diffusion 

coefficient D in the drain boundary is much 
b 

greater than the diffusion coefficient in the 
crystal. Whipple E63 used the Fisher model and 
obtained exact solutions of the average tracer 
concentration profile using Fourier - Laplace 
transforms. He used a constant source condition 
i . e ♦ t 

C (x ? y 1 1 ) =C for t>0 and y =0 

o 

He obtained that the logarithm of the average 

tracer concent ration (In C) varies linearly with 
6/5 

y with a slope that can be related to the grain 
boundary diffusion coefficient by the eouation l 

/ InC V / 4D V f InC \ ' 

V -( -775)1 rX-^r •••«•» 


where } 

n - 

Levine 

boundary 


6 ‘ T5T2 


7 

’ 

and MacCallum L 73 
diffusion eouation 



solved the grain 
directly for the 


average concent ration in polycrystal 1 ine samples 



y 


Free surface 



Grain boundary 


FIGURE 1 - Fisher’s model of a grain boundary. 
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and found that for 


8 >10 the last term in 


Eauation 2.7 is independent of nB~*^ar«d is eaual 
to 0.78. This was later confirmed for Whipple's 
solution by LeClaire E83. 


Suzuoka C9T solved the same problenn but 

assumed an instantaneous source condition ie.» 

C(s<»y»t)=K 6 ( y ) x=0> t=0 

where K is the surface density of the planer 

source. His solution also predicted a linear plot 

6/5 _ 

between In C and y with! — «o.75 

3(nB~ 1/2 ) 6/5 


Thus to obtain the value of D w » experimental 

b 6/5 

values of In C can be plotted against y ,A 
straight line can be fitted to the data points of 
deeper penetrations. The slope of this line can 
be used to determine the value of D w usinsi 


Equation 2.7. 


b 


Zaeschmar C103 treated the grain boundary 
diffusion problem by considering the polycrystal 
as a 3-dimensional brickwork of cylindrical 
grains. He introduced another empirical parameter 
- the filling time of a grain. Oishi and Ichimura 
Clll treated the problem as the homogeneous 
filling of a grain pinned to the grain boundary 



at an arbitrary distance from the surface. As this 
author has used the Oishi-Ichimura geometrical 
configuration for modelling* it will be described 
in detail later. 

2.3 Previous studies of oxygen diffusion in 
alumina 

2.3.1. Lattice diffusivity 

Oishi and Kinsery C123 studied oxygen 

diffusion in single crystal alumina using gaseous 

exchange techniaue. Their plot of log (D> versus 

o 10 

1/T showed a break at 1650 C which was interpreted 

as an extrinsic to intrinsic transition. The high 

temperature data showed that} 

_ , „ ,„-l ,-636 kJ/mole^ 2, , „ B \ 

D ■ 1.9x10 exp( — ) m /sec ...(2.8) 

while the fit to the low temperature data 
suggested that! 

D - 6. 3xl0~ 12 exp(— -- R ^ J ^ mole ) m 2 /sec ...(2.9) 
This was later shown to be an experimental artifact 
sensitive to the surface preparation of the 
specimen. When the samples in the low temperature 
region were preannealed before the diffusion 
anneal* the measured diffusion coefficients fell 
on the extrapolation of the high temperature data. 
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Oishif Ando and Mstsuhiro C133 studied high purity 
vapour grown crystals using the same method as 
Oishi and Kinder*. Their results were in 
approximate agreement with Dishi and Kinder*' s 
results. 

Reed and Weunsch tl43 studied single crystal 
alumina using ion probe. Their data suggest 
thati 

D - 6.4 m 2 / sec . . . (2.10) 

Reddy and Cooper C151 exchanged single 
crystal alumina in enriched 0-18 Sas and measured 
the tracer concentration profile using proton 
activation. Their data fit the eauationl 

D * 2.66x10 ^ ex p ( ■ 6 1 m / S ec ...(2.11) 

Cawley* Halloran and Cooper C163 exchanged 
single crystal alumina in Reddy's gas exchange 
apparatus. The tracer profiles were determined by 
both proton activaton and Secondary Ion Mass 
Spectrometry (SIMS). Cawley found that the 

diffusion coefficient of oxygen was insensitive to 
impurity levels and oxygen partial pressure and 

could be described by} 

-3 .-527 kJ /mole,. 2,_._ ,, N 

D ■ 1.51x10 exp( ) “ ' Bec .-*( 2 . 12 ) 
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Recently Lagerlof * Heuer and Mitchell studied 

the annhilation of dislocation loops by annealing 

deformed cylinders of A1 0 . Based on the 

2 3 

assumption that the self diffusion of oxygen was 

rate controlling they calculated the diffusion 

coefficient of oxygen in undoped alumina asJ 

„ „ . „-4 ,-587 kJ/mole s _2, 

D - 6.8x10 exp( — ) m /sec ...(2.13) 

It should be noted that the data of Oishi and 


Ando* Reddy 

and Cooper* 

Cawley* Halloran 

and 

Cooper* and 

Laserlof and 

Heuer * 

taken by ' 

using 

different 

techniaues * 

agree 

well in 

the 

temperature 

range of 1250 

to 

o 

1800 C. So 

the 

lattice diffusion coefficient 

of oxygen in 

pure 


alumina is well established* 

2*3*2 Grain boundary diffusivity 

Mistier and Coble C173 calculated grain 

boundary diffusion coefficients and widths from 

the results of grain growth* creeps sintering and 

diffusion* Mathematical models of sinterng and 

creep gave D w while grain growth models were 
b 

assumed to give D / w ♦ Thus D and w were 

b b 

calculated separately* The grain boundary width 
for alumina was reported to be between 9-12 nm 


12 



while the dra in boundary diffusion coefficient wss 


reported as! 

n — c /~451 kJ/mole* 2, 

D fc 5.15x10 exp( — ) rn /sec ...(2.14) 

There is a lot of uncertianty in these 
results as the effect of impurities on grain 
drowth was neglected. The model assumed an 

effective diffusion coefficient given by! 

D « f f ■ “b + (1 - £)I> 

where e is the volume fraction of the drain 
boundaries. The above expression does not take 
into consideration the diffusion from the drain 
boundaries into the drains. 

Johnson and Berrin C1B3 used initial stage 
sinterind studies to obtain! 

V ■ 5-^l°~ 4 .xp(-' 53 B !^ / °'° Ie ) » 3 /s.c ...(2.15) 

However their model neglected the contribution of 
surface diffusion to the neck drowth. 

Lessind and Gordon C193 used creep results of 

iron doped alumina samples to obtain! 

„ , „„ , „-7 ,-439 kJ/mole N 3, /0 , 

- 4.38x10 exp( — ) m /sec ...(2.16) 

However they had to guess the rate control 1 ind 
diffusion process. At larde concentrations of 
divalent iron and larde drain size the creep was 
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assumed to be limited ba oxyden drain boundary 


dif fusion. 

Redda and Cooper C15D used proton activation 

on undoped polacrastall ine alumina to obtain? 

. .. « 3 y — 690 kJ /mole.* 2 / i 7 ^ 

D m 4.41x10 exp( ) m /sec ...(2.1// 

b Ki 

The boundary widths assumed were typically 15 to 

40 nm. This assumption is too larde as the drain 

boundary width of alumina is typically closer to 
o* 

10A. There is thus a dreat dedree of disadreement 
between the reported values of drain boundary 
diffusivita and boundary width. 

2.4 Oxidation of superalloas 

2.4.1 Superalloas 

The need to develop materials for hidh 

temperature applications like das turbine endines 

led to the development of superalloas C203. MCrAl 

alloas <M=Fe»Ni»Co) are promisind for use upto 
o 

1200 C. The hidh strendth of these alloys is 

achieved ba solid solution and precipitation 

hardenind while dood oxidation resistance is 

obtained ba the formation of adherent protective 

oxide layers. These alloas owe their oxidation 

# Personal communications with N.L. Peterson of 
Ardonne National Laboratories. 
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resistance to the formation of an alumina scale 

since it is the most stable oxide and has the 

lowest free energy of formation. Alumina is a 

stoichiometric oxide with low diffusion rates and 

is not attacked by either H D or CO . Giddins and 

*? r> 

Pettit 11213 showed that the external alumina scale 
forms for certain compositions of the alloy only. 
The amount of aluminum needed to form a continuous 
scale can be reduced by increasind the chromium 
content as chromium acts as an oxyden Setter. 

2. 4 . 2. Transient scales 

Thoudh the superalloys owe their oxidation 

resistance to the formation of a protective 

alumina scale* it has been found that initially 

transient oxides are formed. , Smialek and 

Gibala C223 studied the structure of transient 

scales formed on pure* Y doped and Zr doped NiCrAl 

alloys. Oxidation of these alloys for 0.1 hour 
o 

at 1100 C produced oriented Y - A1 0 * 

2 3 

a - (Al*Cr) 0 or a spinel like Ni(Cr»Al> 0 
2 3 2 4 

before the random a-Al 0 layer was formed. The 

2 3 

scales were densely populated with internal 
precipitates and there was evidence of plastic 
flow in the underlyind metal due to drowth 
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stresses 


The formation of coherent layers of 


aluminium depleted phases indicated the selective 
removal of aluminium at these short times* 

Golishtlyf Wood and StottE233 studied the 

oxidation of FeCrAl alloys with various Y 

o o 

additions at 1100 C and 1200 C. They observed the 
formation of transient oxides rich in iron and 
chromium in all cases which was followed by the 
steady state a -alumina scale. For the Y free alloy 
it was found that the steady state was reached 
faster if the chromium content or the temperature 
was increased. The amount of transient oxide 

formed was also determined by the specimen surface 
topography since the development of the 
a -alumina layer was less rapid at the base of 
the alloy asperities than at a flat alloy oxide 
interface. 

2.4.3 Mechanism Of Oxidation 

Inert marker experiments by Hindam and 
Smeltzer C243 indicated that alumina scales grow 
inward by oxygen diffusion. Recent 0-18 tracer 
studies by Reddy C251 have also confirmed that the 
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primary transport is oxygen inward and the 
profiles are most consistent with the drain 
boundary transport model* 

GoliShtly et al. C26D observed that the 
growth morphology of alumina scales on FeCrAl 
alloys was different for Y free and Y containing 
alloys. They suggested that in the Y free alloy* 
the scale grew by a combined inward oxygen and 
outward aluminum transport mechanism which caused 
lateral growth within the scale. This was 
responsible for the highly stressed oxide on the Y 
free alloys. The absence of growth stresses and 
lateral growth in scales with Y was attributed 
to the blocking of outward aluminum diffusion. 

Whittle and Stringer C.27J proposed a 
nucleation model to explain the 'Rare earth 
effect 7 ♦ They suggested that the rare earth 
additions acted as nucleation sites for all the 
oxide phases* thus producing a finer grained 
initial oxide. This caused the steady state 
oxidation to be reached more rapidly and at a 
lower aluminum or chromium content. 
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2.4.4 Oxide adherence 


As discussed earlier* the introduction of 
small amounts of rare earth elements to the 
superalloys caused a dramatic improvement in the 
scale adherence. Some of the mechanisms put 
forward to explain this phenomenon are as 
follows £27*281] i i [ 

a) Y (or other rare earth elements) improve the 
chemical bond at the metal-oxide interface as they 
have a better affinity for oxygen than the 
alloying elements. 

b) Yttride precipitates or yttria particles can 
induce mechanical pegging of the oxide scale on 
the substrate. 

c) Y promotes the formation of an inner yttria 
layer which acts as a diffusion barrier and also 
decreases the mechanical and thermal stresses 
between the oxide and the alloy. 

d) Soluble Y in the alloy can trap the metal 
vacancies and conseouently avoid vacancy 
coalescence at the metal oxide interface. 
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e) Y present in the oxide scale can infuence oxide 


defect mobility and conseauently the diffusion 
rate in the oxide scale * oxide thickness and 
plasticity. 

Delaunay et al. studied the mechanical 
stresses generated in the oxide scale during the 
oxidation of superalloys. They concluded that Y 
improves the oxide adherence by trapping the 
vacancies and preventing them from precipitating at 
the metal oxide interface E29D.lt may also improve 
the oxide plasticity by increasing the 

concent rati on of oxygen vacancies E30D. 

2.4.5 Growth kinetics and microstructure of 
FeCrAl-Zr alloy 

Smialek studied the microstructure of the 
oxide grown on the FeCrAl-Zr alloy used by this 
author. He saw monoclinic zirconia precipitates 
at the grain boundaries at the gas surface. There 
were deep depressions in the large grains at the 
metal-oxide interface which represented areas 

where the grain boundary regions had advanced 
faster than the interior of the grains due to 
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o 

short circuit diffusion C31f32»33D. At 1100 C» 
the growth rste of the oxide was found to obey the 
law? 

0.46 

X = (K t) ...(2.18) 

p 

Where X is the oxide thickness and t is the 
oxidation time. The grain growth in the alumina 
scale was found to obey the law; 

0.2 

Grain size = (K t) ...(2.19) 

g 

There was a variation of grain size and shape 
across the scale. The grains near the gas surface 
were spherical and small in size while the grains 
near the oxide-metal interface were larger in size 
and columnar in shape. 

2.5. Previous 0-18 sudies 

0-18 has been successfully used as a tracer 
to identify the growth mechanisms of oxide scales. 
The double oxidation techniaue consists of 
successive oxidation of the alloy in natural 
oxygen and 0-18 tracer. The concentration profile 
of the tracer is determined by some convenient 
techniaue like proton activation or SIHS. The 
tracer can then be used to identify the regions of 
new growth. The shape of the tracer concentration 
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profile can provide insidht into the mechanism of 


oxide slrowth. 

0-18 has been succesfully used as a tracer by 

many workers. Lees and Calvert E341 used the 

double oxidation techniaue on FeCr alloys in 1976. 

In 1978? Sheasby and BrownE351 used SIMS to obtain 

0-18 concentration profiles in CoO scales. In 

1980» Hartley et al. 1261 at Harwell oxidized 

18 

FeCr alloys in a C 0 atmosphere. They used 

2 

proton activation to obtain the tracer 
concentration profiles. In 1982 Lees et al. 
studied the mechanism of corrosion of low silicon 
E37D and hidh silicon E381 ferritic steels in hiah 
pressure CO usinsi 0-18 tracer. 

ReddaC 153 used the double oxidation technioue 
to measure oxyden dif fusivities in various oxides. 
He used the proton activation technioue to measure 
the tracer concentration profiles. Reddy also 
studied the oxidation mechanism of NiCrAl alloys 
by the double oxidation techniaue. Very recently 
Cawley studied the self diffusion coefficient of 
oxyden in alumina usind 0-18 tracer. He 
determined the tracer profiles by both proton 
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activation and SIMS 


Cawley also studied the 


thermal oxidation of silica and silicon carbide 
using the same technioue. He ^as able to obtain 
good fits of his data to the theoratical 
calculations if he assumed different diffusion 
coefficients for his near surface and deep 
profiles « 

2.6. SIMS 

In SIMS C393*the sample is bombarded with 

+ 

high energy primary ions (Ar in this case). These 
primary ions penetrate into the solid where binary 
collisions with target atoms occur. The primary 
atoms come to rest at the so called penetration 
depth and are thus implanted in the solid. These 
primary ions cause the target atoms to be 
displaced from the surface due to various 
generation target recoils. These target or 
secondary ions are passed through a mass analyser 
where they are separated by mass to charge ratios. 
An electron beam is also focussed on the sample to 
keep it neutral. 
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The number of sputtered secondary ions of a 
given type is a complex function of the surface 
concentration* electronic structure and sputtering 
conditions* However* ss all the variables for 
isotopes are almost the same* the ratio of counts 
of isotopic ions gives the true value of isotopic 
concent rat ion « 

The resolution of SIMS is limited by 
non-uniform sputter etching which may yield 
sputtered species from a range of depths C40.411, 
This problem is particularly severe at the edge of 
the crater. To avoid this problem the raster 
gating technique is used where the beam scans over 
a large area and detects secondary ions from the 
central portion of the crater only. Another 
problem faced when sputtering insulators is that 
of charging due to ion removal. If the charging 
is stable it can attenuate the secondary ion 
signal C42D. The electron beam can neutralize the 
secondary ions* causing further complications. 
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The advantage of using SIMS over proton 


activation is its depth resolution which is less 
o 

than 10A. Also there is no limit on the depth of 
the sputtered crater exept for the machine time* 
SIMS also provides the data in a very convenient 
form and can monitor upto 9 elements at a time* 

However? the problems of non uniform 
sputtering? especialla in polacrastalline samples 
can severela limit the resolution of SIMS and can 
actualla make the data worthless. This point will 
be discussed in detail later. The conversion of 
counts to concentration with the help of standards 
can be extremely unreliable as the ion yield 
depends on? among other things? the sample surface 
preparation? the incident angles of the ions and 
the sample microstructure. So it is extremely 
important to apply stringent tests on the Quality 
of data before any Quantitative information is 
obtained from it. 
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CHAPTER 3 


THEORY 

3.1 Statement of the problem 

0-18 has been successfully used as a tracer to 
identify the growth mechanisms in oxide scales. 
The double oxidation techniaue is used to identify 
the diffusion processes involved. In this 
techniaue the alloy is seauentially oxidized* first 
in natural oxygen and then in 0-18 tracer. The 
oxide grown during the natural oxygen oxidation is 
referred to as the old oxide while the oxide grown 
during the tracer oxidation is referred to as the 
new oxide. The position of the new oxide growth is 
determined from the concentration profile of the 
tracer across the oxide scale. The tracer 
concentration profile is obtained by SIMS. 

It has been shown that the primary transport 
mechanism operative during the growth of a -alumina 
scales on superalloys is inward grain boundary 
diffusion of oxygen C15*321 .Thus the protective 
oxide scale has been found to be primarily inward 
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arowina 


However » as the oxygen tracer diffuses 


down the grain boundaries » it exchanges with 0-16 
existing in the grains. To get any Quantitative 
information from the tracer concentration profile * 
it is necessary to model this process of inward 
tracer grain boundary diffusion with exchange. The 
advantage of a Quantitative model is that the grain 
boundary and the lattice diffusion coefficients can 
be obtained from the shape of the tracer 
concent ration profile. The effect of structural 
parameters like the oxide grain size on the tracer 
profile can also be studied. It can also be used to 
determine if the scale is inward or outward 
growing. In the case where scale growth is a 
combination of both the mechanisms* a simple mass 
balance test can be used to separate the two 
processes. To model the diffusion process* Dishi 
and Ichimura's geometrical configuraton of the 
grain boundaries and grains have been used. 

In the model* the shape of the tracer 
concentration profile depends on the three 
parameters? 

1) the grain size r* 2) the grain boundary 
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diffusion 


coefficient D and 3) the lattice 
b 

diffusion coefficient D. By choosing the 
appropriate values of the three parameters * a 
family of curves can be generated which define a 
range of 0-18 profiles which are auanti tatively 
compatible with the inward oxygen diffusion 
mechanism. Conversely* given an experimental 
tracer concentration profile* it should be possible 
to pin down the values of the three parameters. It 
should be noted that the grain size of the oxide 
scale can be obtained by SEM studies. A good 
estimate of the grain boundary diffusion 
coefficient can be obtained from the growth 
Kinetics of the oxide while values of the lattice 
diffusion coefficient can be obtained from the 
literature. Thus there is independent information 
existing for all the three model parameters. A 
major assumpton made during the modelling of the 
diffusion process is that the oxide scale grows in 
a planar fashion. 

3.2 Expected Tracer Profiles 
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Figure 2 shows the schematic 0-18 

concentration profiles for various transport 

mechanisms. If the outward aluminum transport is 

the primary mechanism the scale will be outward 

growing. The tracer concentraton profile will be a 

18 

step function due to the formation of A1 0 at the 

2 3 

gas-oxide interface. If the transport mechanism is 

inward oxygen grain boundary transport* the scale 

will be inward growing and the concentration 

profile will be a step function due to the formaton 
18 

of A1 0 at the metal-oxide interface. This is 
2 3 

however an ideal case. When the tracer in the 
grain boundaries comes in contact with 0-16 in the 
grains of the old scale* exchange will take place 
and the concentration profile will look like Figure 
2c). Finally if the diffusion mechanism is a 
combination of outward aluminium and inward oxygen 
transport the tracer profle would look like Figure 
2d) . 

3.3 Mass Balance 

Given a tracer concentration profile across 
the oxide scale a simple mass balance test can be 
used to determine the growth mechanism. As shown 
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ORIGINAL PASS ® 
Qf* POOR QUALITY 



OUTWARD A1 TRANSPORT 



WWARD OXYGEN GLB. TRANSPORT 
WITHOUT EXCHANGE 


a) 


b) 



wrm EXCHANGE 



INWARD OXYGEN TRANSPORT 


c) 


d) 


FIGURE 2 - Schematic tracer concentration profiles 
for various transport mechanisms. 
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in Figure 3? a normalized concentration profile can 


be plotted such that the X-axis is a plot of the 
fraction of the total scale thickness while the 
Y-axis is the plot of the concentration of the 
tracer as a fraction of the maximum tracer 
concentration in the scale. 

We know that evera atom of 0-18 in the 

oxide scale sot there when the sample was exposed 

to the tracer. Thus the number of lattice sites 

occupied ba the tracer atoms in the oxide scale 

must be the same as the number of oxasen lattice 

sites created in the oxide during the tracer 

oxidation. This concept can be used to calculate 

the thickness of the new oxide Y as shown in 

new 

Figure 3b). Y is so chosen? that area A in 
new 

Figure 3a) (amount of 0-18 in the total oxide) is 
eoual to the area B in Figure 3b) (amount of oxasen 
in the new oxide). Having established the new 
oxide thickness? it can be determined if the growth 
mechanism is primarila inward oxasen transport or a 
combinaton of outward aluminum and inward oxasen 
transport. If the diffusion mechanism is inward 
transport of oxasen through the grain boundaries - 
with exchange? the scale is inward growing. Then 
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a) 


|“ Y old ^newf— 



c) 


FIGURE 3 - Determination 
mass balance. 



d) 


of transport mechanisms by 
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the number of 0-18 atoms retained in the old oxide 


must be the same as the number of 0-16 atoms found 

in the new oxide as they are both involved in the 

same exchange process* Thus if area C is eoual to 

the area D in Figure 3c) the transport mechanism is 

inward oxygen diffusion* However* if area C is 

greater than area D some of the tracer atoms in C 

did not get there due to exchange but rather due to 

outward aluminium diffusion leading to the 
18 

formation of A1 0 at the gas-oxide interface* In 
2 3 

this case* it is possible to determine* at least in 

principle* what fraction of the scale growth is due 

to each mechanism. Thus Y can be split up into 
in out new 

Y and Y such that area E=area F in Figure 3d), 
new new 

This allows the growth kinetics to be split up into 
inward growth and outward growth and both these can 
be separately analysed. So with profiles of 

sufficient Quality* one can obtain detailed 
information by simple mass balance analysis* It 

should be noted that this is model free. 

It should be mentioned that the mass balance 
criterion is valid if the number of tracer atoms in 
the grains at the oxide-gas interface due to direct 
surface exchange with the gas is neglected. This 
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is s valid assumption because the drain boundary 


diffusivita is much higher than the lattice 
diffusivita in the oxide arid the penetration depth 
of the tracer atoms due to volume diffusion at the 
surface is negligible. 

Much more informaton can be obtained from the 
actual shape of the tracer concentration profiler 
if the actual eauations describing tracer 
incorporation during scale growth could be directly 
solved. Even in the simple case of inward oxygen 
diffusion it is a formidable problem of combined 
lattice and grain boundary diffusion with a moving 
boundary under an imposed field. However a 
relatively simple but still realistic model permits 
the diffusion eauations to be solved. 

3.4 Geometrical Configuration 

Figure 4 shows an idealised drawing of the old 
scale with hexagonal grains arranged in a close 
packed two dimensional network. These grains are 
formed during the 0-16 oxidation and contain 
natural oxygen atoms. When the sample is exposed 
to the tracer* 0-18 diffuses rapidly down the grain 
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ORIGINAL PAGE ISf 
OF POOR QUALITY 


0-18 GAS 



FIGURE 4 - Schematic of grain boundary diffusion of 
tracer with exchange. 


34 


boundaries to form new oxide at the oxide-metal 


interface* However* while diffusing down the grain 
boundaries the tracer comes into contact with the 
grains which have 0-16 in them. This causes a 
tracer exchange between the grain boundaries and 
grains. It should be noted that the grains close 
to the gas surface have more tracer in them due to 
greater exchange time than the grains at the 
old-oxide new-oxide interface. 

To model this process of grain boundary 
diffusion of tracer with exchange* the 
Oishi-Ichimura geometrical configuration has been 
chosen as shown in Figure 5. The eolycrustal 1 ine 
oxide is assumed to consist of eaual sized 
spherical grains and to possess a grain boundary 
diffusion coefficient that is much greater than the 
lattice diffusion coefficient. Diffusion of the 
tracer proceeds from the gas surface preferentially 
along the grain boundaries and from the grain 
boundaries radially inward towards the center of 
the grains. The grain boundaries of thickness 
are assumed to be perpendicular to the oxide 
surface « The grains of radius r are assumed to be 
attached to the grain boundaries at any arbitrary 
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0-18TRACER 


OLD OXIDE 

(FORMED DURING 
0-18 OXIDATION) 



OXIDE GRAIN 


NEW OXIDE 


iH 


GRAIN BOUNDARY 
TOTsJy+Jy* 


OXYGEN USED TO 
FORM NEW OXIDE 


METAL 

FIGURE 5 - Geometrical configuration of oxide grain and grain boundary. 



depth in the oxide and are assumed to be completely 
surrounded by the grain boundaries* Basically* the 
grains act as sinks for the tracer and are attached 
to the grain boundary diffusion path. The 
complicated two dimensional combined lattice and 
the boundary diffusion problem is transformed into 
a simple one dimensional diffusion problem in the 
grain boundary • with the grains acting as 

distributed sinks. The strength of the sinks is 
given by the well known case of diffusion into a 
sphere . 

It should be noted that when the oxygen flux 

TOT 

J enters the grain boundary at the gas surface* 

it is all 0-18. However* as it comes in contact 

18 

with a grain* there is a flux of tracer* J » 

* * 

A 

into the grain due to exchange. However as the 

number of oxygen sites in the grains are fixed* 

16 

there is an eaual and opposite flux J of 0-16 

* \ 

atoms out of the grain and into the grain boundary. 
Thus as the tracer diffuses down the grain boundary* 
it becomes more and more diluted in 0-16 due to 
exchange. The grains act as a sink drawing 0-18 
atoms from the grain boundary and as a source 
supplying 0-16 to the grain boundary. The grain 
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oxide grain 


grain boundary 


new oxide 


FIGURE 6 - Coordinate system for the model. 
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boundary oxygen flux at the old oxide-new oxide 

'Is 

interface which is used to form new oxide is 
actually a mixture of tracer and natural oxygen. 
This suggests that the isotopic content of the new 
oxide drains preserves the information on the 
amount of exchange occuring during the transport of 
the tracer through the old oxide. 


Figure 6 shows the coordinate system used in 

the mathematical treatment of the diffusion 

process. The list of symbols used are as follows. 
18 2 
J t Flux of tracer down the grain boundary (kg/m 
y 

sec ) 

18 

J { Flux of tracer out of the grain boundary (kg/m 
x 

sec) 


D i Grain boundary diffusion coefficient (m /sec) 
b 

n 

JD ♦ Lattice diffusion coefficient <m /sec) 


R i Gas constant (Joules/K mole) 


T I Temperature (K) 
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y° I Oxidation chemical potential < Joules/K mole) 


Y * Oxide thickness (m> 


K J Parabolic rate constant for oxide growth(m 

p 

/sec ) 


t X 0-18 oxidation time (sec) 


x X 0-16 oxidation time (sec) 


d> J Tracer concetration in the grain boundary 

y 3 

at depth y ( ks/m ) 

3 

C X Tracer concent ration in the grain( k d/ m ) 


C 1 Initial tracer concentration in the grain 
i 3 

( kg/m ) 

3 

C ♦ Average tracer concentration at depth y(ks/m ) 


tf X Grain boundary width (m) 


r J Grain size (m) 


CJ Average tracer concentration in oxide 
g 3 

grain(kg/m ) 
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3.5 Mathematical Modelling Of Diffusion 


The flux of tracer down the grain boundary 


18 


(J ) is due to two gradients* the oxidation 
a . O' 

chemical potential sradient( vy°-^- > which 

drives oxidation and its own concentration 

30 

sfradient( y ) due to isotopxc mxxins. 

9y 

D. 0 A 30 

J y “"RT^ "T~ " D b3T ...(3.1) 

The diffusion coefficient D is different from the 


b 

self diffusion coefficient of oxygen as Jumps in 
the direction of the gradient are favoured. There 
is also a difference in correlation coefficients. 
However this difference in diffusion coefficient is 
small C433. 


The flux of tracer out of the grain boundary 

18 

and into the grain (J ) is proportional to the 

x 

concent ration Gradient of the tracer in the srain 
surface and the lattice diffuson coefficient D. 


-18 

J x A ...(3.2) 

The srain boundary is assumed to be thin enough to 


neglect the variation of in the x direction. 

Fick's second law gives the time dependence of the 


grain boundary tracer content. 

18 



where & is a sink, term representing the rate of 
tracer loss into the grain* g is given by E53. 


6 - ! 4 s 

Substitution of the Eauation 3.2 
gives* 


B - - 


2D . 8C. 
6 1 dx } 


...(3.4) 

in Eauation 3.4 
...(3.5) 


Thus the time dependence of the tracer content in 
the grain boundary is given by substituting 
Eauation 3.1 and Eauation 3.5 in Eauation 3.3. 


3 $ 

at 


V u ° V „ iV 2D, 3C. ,, .. 

rtT" bT b^2 + ) 


The last term of the Eauation 3.6 can be evaluated 
if the concentration gradient of the tracer at the 
grain surface is known. The increasing rate of the 
grain boundary concentration is much larger than 
the increasing rate of the tracer concentration 
within the grain for a given time* since the grain 
boundary diffusion coefficient is much larger than 
the lattice diffusion coef f icient . Thus the 
concentration variation within the grain can be 
considered for a constant 4> y with a relatively 
small error. Thus the problem of tracer exchange 
reduces to the standard problem of diffusion into a 
sphere having a constant surface concentration <f y ♦ 
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kv&nal fwa a 
OF. POOR QUALO'V: 

The standard solution to this eouation is siven in 

Eauation 3.7. 

C-C. 


1 + 


2 r 


♦y^l ’ (r - x) n-l ° 


I<rl) n 6ln Elilrii ^(ZBnflA, ••■( 3 . 


-) 


7) 


Differentiation of Eauation 3.7 by x slives Eauation 
3.8f an expression for the tracer concentration 
Gradient at the arain surface (if 5> « r). 


3c 00 2 2 

( 3x^ 6 “ -<K- C i>T 1 expt^nlL.t) 

n*l r 


...( 3 . 8 ) 


Substitition of Eauation 3.8 in Eauation 3.6 sives 
Eauation 3.9. 


3$ B4> 

y b 


3 2 <p 


AD. 


2 2 


St ~ RT Y 9 y + ^TT 21 ” 1 ex P (~ Dn 2 t ) -..( 3 . 9 ) 


3 y 


n-1 


This is a differential eauation involving only the 
strain boundary concentration and its variation in 
the y direction and time. The scale thickness is 
time dependent. So if a parabolic drouth rate is 


assumed with K as the parabolic rate constant? 

p 

1/2 

Y - (K p (t+T)) '*■ ...(3.10) 

where t and T are the 0-18 and 0-16 oxidation times 


respectively. Eauation 3.9 can thus be written as.' 

o 


3ct> 

— * - n 


sV 




9<f> 


2 2 
-Dn tt t. 


- $ _r 


3t b 9 y 2 RT(R p (t+T )) 1/2 9y ‘ w n -1 

4DC , ® ~2 2 

+ Z exp( 2 — ) ...(3.11) 

This differential eauation can *be easily solved 


43 



numerically to give the grain boundary tracer 

concentration 4>y(y»t). Note that depends upon 

three transport parameters D t 6 and I)r the oxide 

b 

microstructure through the grain size r and the 

oxide growth kinetics through the parabolic 

constant K ♦ When Eouation 3.11 is divided by the 
p 3 

total oxygen concentration in kg/m t which is the 

same throughout the oxide scaler the tracer 

3 

concentration can be converted from kg/m to 
fractional oxygen concentration. Thus Cr C and 4> y 
all become unitless and are actallyJ 

kg of tracer per cubic meter 
kg of oxygen per cubic meter 

As the SIMS determines the tracer 
concentration profile by the sectioning techniauer 
the average concentration of the tracer at any 
depth y should be determined. The average tracer 
concentration in the oxide grain at depth y at a 
time t is given by CUD 

V t ( W + C i ...(3.12) 

where M and are the total amounts of tracer 
t 

diffusing into the grain at times t arid » r 

M. 

respectively. _E is not a function of the grain 

00 
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boundary diffusion coefficient D but only of the 

b 

lattice diffusion coefficient D and is given by the 


standard form till} 
M 


CO O O 

t , 6^1 ,-Dn TT t x 

5 T " 1 2 r , — r~) 

00 Ti n**l n r 


..(3.13) 


Thus C can be written as? 

a 

— 00 

c 


g 


l ~ exp( ^ P ))(<fr -C ) + C j 
v n»l n r y 


. . (3. 1A) 


If the contribution of the tracer in the grain 
boundary is considered* the average tracer 


concentration at depth y and time t is given byJ 

4 tt r 3 C + 4 tt x 2 | 4> 

c (y,t) - 8 L-Z ...( 3 . 15 ) 


f ir r 3 + 4, r 2 | 

Due to the rasterinS technioue C(y*t) is the 
auantity measured by SIMS! 


3.6 Boundary Conditions 

Eauation 3.11 is solved by a numerical 
solution using the finite difference method. 

<y*t> is determined by using the first central 
difference form of the deriviatives £453. If X is 
the thickness of the old oxide and Y(t) is the 
total thickness of the oxide at any time t»the 


initial conditions are? 


4> - C i , y<X , t-0 
C ■ , 0<x<r , y , t*0 


..(3.16) 
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It is assumed that the natural 


isotopic 


concentration of the tracer is negligible i.e.» 

— 0 • • • (3.17) 

The boundary condition at the oxide-gas interface 

1S * * y “ C 6 » t>0 » y*0 ...(3.18) 

where C is the tracer concentraton in the gas 
s 

during the 0-18 oxidation. C is assumed to be 

s 

constant during the diffusion. 

The oxide-metal interface is a moving 
boundary. As new oxide is being formed at the 

oxide-metal interface due to the flux of tracer 
through the grain boundaries* the tracer 

concentration in the new oxide grains being formed 
is the same as the grain boundary tracer 

concentration at the metal-oxide interface. Figure 
7 shows the new oxide growth Ay during the time At. 

The boundary conditions at this interface are 5 

C(0<x<r, Y<y<Y+AY, t-t+At) - 4> y (y* Y » t-t+At) 
4> y (Y<y<Y+AY, t-t+At) - 4> y (y-Y, t-t+At) ...(3.19) 

where limit » AY, At -*• 0 

In the following section some of the auantities 
necessary for the numerical solution of Eauation 
3.11 will be determined. 
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FIGURE 7 - Boundary conditions at the oxide-metal moving boundary. 



3.7 Determination of oxidation chemical potential 
gradient 

Let the gas oxide interface be denoted as 1 

and the oxide metal interface be denoted as 2 as 

I 2 

shown in Figure 8« Let and be the oxygen 

u 2 0 2 

chemical potential in the oxide at 1 arid 2 
respectively. At 1 the chemical potential of 
oxygen in the gas must be the same as the chemical 
potential in the oxide C463 i.e.» 


V, 


,gas 


"2 2 

As the sample is exposed to a closed atmosphere of 
pure oxygen during the tracer oxidation » 

P? a8 - 1 


This gives 


1 


- P® 38 - RT la P 8 as - 0 ...(3.20) 


2 ~2 v 2 

i.e«» we choose the standard state as 1 atmosphere 

oxygen gas. At 2 the oxygen chemical potential 

must be the same at the oxide and the metal i.e.» 

2 metal 

V - V 


At 2 the reaction! 
2A1 


(s) + 2 °2(g) 


A1 2°3(b) 


...(3.21) 


A ft 
a o 




takes place which has a free energy change of 


An eauilibrium constant K for this reaction can be 
defined as! 

AC° . -RT In K - -RT la ^ n -(3.22, 

(a Al> (a 0 2 > 

“at n can be assumed to be unity. Thus Eouation 
2°3 

3*22 becomes! 


£ etal « exp(-ll^i)(a )- A/3 
0 pv 3 RT ^ a A1 '» 


Since 
dives ! 


2 

metal 


...(3.23) 


- RT In a 


metal 


r Eouation 3.23 


.metal _ 2 i 4 „oi 4 


~ T l AG °l “ T RT In (a..) 


...(3.24) 


'2 - A1 

In the absence of any thermodynamic activity data* 

the metal can be assumed to be a sindle phase ideal 

solution. Then a is the atomic percent of 

A1 

aluminum in the alloy f EA13. Thus the Eouation 
3.24 becomes! 


metal 


- - 3 |AC | - | RT !n[Al3 ...(3.25) 


2 ”2 

Thus the oxidation chemical potential difference 
across the oxide scale is given by! 

0 “ \ “ S' f |AGW| + 3 
Substitution of values of and CA13 in Eouation 
3.26 dives! Au° - 787915 J/mole at 1100°C. 


4 / 


P 0_ " p 0 “ 3 l AG °i + 4 RT ln C A1 ^ ...(3. 


26) 
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3.8 Relationship between D and K 

b * 

If the mechanism of oxide growth is primarily 
inward oxygen diffusion through the grain 
boundaries? the grain boundary diffusivity and the 
growth rate kinetics can be related by considering 
a simple model. 


The sample is exposed to natural oxygen and an 

oxide scale is grown. The oxide is assumed to have 

cubic grains of dimension r surrounded by the grain 

boundaries of thickness 5 as shown in Figure 9. 

The cross section of the sample itself is assumed 

to be a souare with each side being n grains long. 

Let the thickness of the oxide at any time t be Y 

and let the scale grow by an amount AY in time dt. 

o 2 

Then the flux of oxygen J <kS/m sec) down the grain 
boundaries is due to the oxygen potential gradient. 


The amount of oxygen transferred in time 

o 

through the grain boundaries? U (kg) is given by! 

T 

W° “ J°. area of grain boundaries . time 
the shaded area in Figure 9 is the grain boundary 


area through which the oxygen is transported (6 «r ) . 

Dl C A. O 


b Au „ 2 r 
W T “ RT~ Y~ • 2n r $ • At 


...(3.27) 
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FIGURE 9 - 



nfr+6) 


grain boundary 



oxide 

grain 


Cross section at metal-oxide interface 11* 


Schematic of model geometry used to relate scale growth kinetics to grain 
boundary diffusion. 



as o 


This amount of oxygen transported must be the 

same as the amount of oxygen incorporated in the 

o 

net growth of oxide* W <kg). 

N 

■ volume of new growth in time At . concentration of 0 in scale 


V° - (nr) 2 . AY . C 


. . ( 3 . 28 ) 


where C is the oxygen concentration (k-s/m ) in the 
o o 

scale. As W » U » Eauatiori 3.27 and Eauation 
T N 

3.28 dives? 

D, C Ay « 2 2 

— 2n r 6 At - n r AY C 
RT Y 


26D,Ay° AY 


.( 3 . 29 ) 


AY 

At 


RT r Y At 

is the growth rate of the oxide. If a 
parabolic growth rate is assumed with growth rate 
constant K_ 


Y - (K p t) 


1/2 


- 1/2 

dt 2 p • 

Combining Eauations 3.29 and 3.20 gives? 


( 3 . 30 ) 


®b ' 75 

Ay 


..( 3 . 31 ) 


Substitution of appropriate values in Eauation 

o -16 2 

3.31? r=0.3 microns* T=1100 C * K =1.2x10 m /sec 

o p 

* 5 cio A * Ay° =78915 J/mole gives? 
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...(3.32) 


11^= 1.3x10 m /sec 

This gives a good estimate of the value of D b to 
be used in the numerical solution of the Eouation 
3.11. 

3.9 Tracer concentration profile in the new oxide 

The tracer concentration profile in the new 
oxide is obtained by the principle of mass 
balance. It is assumed that the amount of exchange 
of 0-16 in the oxide drains at the das surface with 
the 0-18 in the tracer due to volume diffusion at 
the das-oxide interface is negligible. Then for an 
inward growing scale once the sample is exposed to 
pure tracer* the amount of 0-16 in the entire oxide 
must remain constant. This means that the amount 
of 0-18 retained in the old oxide must be the same 
as the amount of 0-16 ,in the new oxide. Let X be 
the thickness of the oxide after an 0-16 oxidation 
for time t and let Y be the oxide thickness after a 
further tracer oxidation for a time t. Let an 
additional tracer oxidation for a time dt cause a 
new oxide growth of dY as shown in Figure 10. The 
tracer concentration profile till depth Y at time 
t+dt is also shown. 
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Ares A denotes the amount of tracer in the old 
oxide and ares B denotes the amount of 0-16 in the 
new oxide grown till time t. Let A-B=E» Then by 
mass balance* the amount of 0-16 contained in the 
new growth dY must be E. As mentioned in Section 
3.6» the tracer concentration in the new growth is? 

♦y ■ (yY. t-t+6t) 

which gives the 0-16 concentration in the new 
growth as: (i _ yy. Y> t «t+At) 

Thus the growth dY during time dt can be calculated 
as J 

dY = (l-$ y (y=Y^t=t+At)) •••<3-33) 

3.10 Flow chart of computer solution 

Eouation 3.11 is solved on the computer using 

the finite difference techniaue. The old oxide is 

divided into 19 eaual parts with the oxide-gas 

interface being called 1 and the metal-oxide 

interface being called 20» as shown in Figure 11, 

Each part is considered to be 1 unit. The growth 

rate of the oxide is determined by mass balance as 

discussed in Section 3.9, The value of D is 

b 

chosen such that the growth rate is compatable 
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with the parabolic growth rate kinetics* The flow 
chart for the computer program is shown in Figure 
12* The computer program is listed in appendix A. 

3*11 Computer plots of tracer profiles 

Figures 13-15 show the theoretical tracer 

concentration profiles in the old oxide for various 

values of lattice diffusivita* grain boundary 

diffusivita and grain size. A higher value of 

lattice diffusivita and causes a fast exchange rate 

to be set up. The grains in the old oxide are filled 

up auickly. A higher value of the grain boundary 

diffusivita causes more tracer %o be retained in 

the grains of the old oxide as the tracer diffuses 

auickly down the grain boundaries increasing the 

exchange time for the grains. A smaller grain size 

increases the surface area for exchange and causes 

more tracer to be retained in the old oxide. The 

variation of the tracer profiles when each 

parameter is varied is different. Thus 

theoretical la given an experimental tracer profile* 

it should be possible to come up with uriiaue values 

of the three parameters to fit the curve. It 

should be noted that the values of D and r are 

b 
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INITIALISE 

/ y (T=0)=0: N=l,21 

C(T=0) =0: N=l,20 
NCOUNT=20 

E 


{ T=T+TIN’C 


| CALCULATE / y (T) : N=1 .NCOUNT 


yes 


CALCULATE C 

: N=l, NCOUNT 



| PRINT C : N=1 .NCOUNT 



yes 


HAS 

["THE OXIDE GROWN 
BY 1 UNIT ACCO 
RDING TO PARA- 
BOLIC KINETIC 


CALCULATE SCALE GROWTH 
(IT SHOULD BE COMPATIBLE 
TO PARABOLIC KINETICS) 



CALCULATE C:N=1 .NCOUNT 
PRINT C: N=1 .NCOUNT 


I 


0 y (NCOUNT+1) (NCOUNT) 
C (NCOUNT+1) =/ y (NCOUNT) 


INCOUNT « NCOUNT* 1 


CALCULATE C : 
N=1 , NCOUNT 
ZZ I 

PRINT C : 

N=1 , NCOUNT 

~ ■ z 

CALCULATE 
OXIDE GROWTH 


no 



FIGURE 12 - Flow chart for the computer program. 




0-16 oxidation time- 18 hours 
0-18 oxidation time- 14 hours 
at 1100°C ,, - 
D b =1.3xl0'- LO m z /sec 
grain size=0.3 microns 



FIGURE 13 - Effect of lattice diffusion coefficient 
on the tracer profile in the old scale. 



F; POOR QUAU7 


0-16 oxidation time-18 hours 
0-18 oxidation time-14 hours 
at 1100°C 7n 7 
D=5.0xl0‘ mVsec 
grain size=0 . 2microns 



FIGURE 14 - Effect of grain boundary diffusion coeffi- 
cient on the tracer profile in the old 
scale . 
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not independent » but are related to the K value by 


Eauation 3.31. Thus if one of the above parameters 
is fixed? there is a uniaue value of the other 
parameter to satisfy the parabolic growth rate 
kinetics. Figures 13-15 are plotted to show how 
sensitive the profiles are to the three parameters. 

Figure 16 shows the computer simulation of how 
the tracer concentration profile develops accross 
the entire oxide scale when the scale grown by a 
18 hour 0-16 oxidation followed by a 1 4 hour 0-18 
oxidation. The values of the three parameters 
chosen are! 

-16 2 -20 2 
D = 1.41x10 m/sec> D = 10 m/sec 
b 

r = 0.3 microns 

The value of the grain boundary width chosen is 
o 

10A * 


Here we can see the versatility of the model. 
The oxide microst ructure can be incorporated in the 
program by feeding in the real grain size variation 
across the scale. In case the lattice diffusion 
coefficient varies across the side* that too can be 
incorporated in the program. The actual value of 
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the oxide growth Kinetics can be experimentally 

found and the real value of Y(t) can be used in the 

program. The value of D can be so adjusted that 

b 

the growth rate predicted by mass balance is the 
same as the experimental growth rate. So long as 
the oxide is inward growing? the tracer profile for 
any growth kinetics and microstructure can be 
modelled . 

It should be noted that the model is based on 
a Planar oxide growth mechanism and any 
Quantitative information derived from the 
experimental profiles assumes good resolution and 
high Quality of the data itself. In case the oxide 
metal interface is irregular? it is difficult to 
identify the interface for the purpose of mass 
balance. The same problem arises if the sputtering 
itself is non uniform and information at the 
metal-oxide interface is being gathered 

simultaneously from both metal and oxide grains. 
Non uniform sputtering can severly limit the 
resolution of the data such that no Quantitative 
information can be obtained from the shape of the 
profile. 
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CHAPTER 4 


EXPERIMENTAL PROCEDURE 

4.1 Sample Preparation 

The FeCrAl-Zr alios used in the experiment 

# 

was obtained through NASA Lewis Research Center. 
The alios used was Fe-21.9Cr-5.8Al-0.22Zr and was 
obtained as cold rolled strips about 1.5mm thick. 
Samples of size 6mmx6mmxl . 5mm were cut out bs a 
Buehler diamond saw. The specimens were then 
mounted and polished on 240»320 and 600 grit 
polishing paper followed bs 15»3 and 1 micron 
diamond paste polishing. Some of the samples were 
prepared at NASA and were polished bs 600 grit 
polisins paper followed bs 3 micron and 1 micron 
diamond paste polish. . 

4.2 Exchange Apparatus 

Figure 17 shows a schematic of the apparatus 
used for the gas exchange process. The tracer is 
kept in a zeolite storage tank. 

* Courtses J.Smialek. 
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FIGURE 17 - Schematic of gas exchange apparatus. 
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The sample is placed in a Quarts tube with 


the polished side up. Oxygen flow into the Quarts 


tube can 

be 

controlled 

by valve 

5 and 

the 

stopcock . 

The thermocouple vacuum 

gauge measures 

the pressure 

inside the 

system 

when it 

is 

evacuated. 


The pressure gauge 

measures 

the 

pressure 

in 

the system 

when the 

tracer 

is 


introduced. Normally the tracer pressure in the 

system is maintained slightly above the 

atmospheric pressure such that no air leaks into 
the system during the tracer anneal. The sample 
is annealed by placing the Quarts tube in a single 
zone horizontal Lindberg furnace. The temperature 
of the furnace is monitored by a F't~6%/Rh/F‘t-30%fth 
thermocouple . 

4.3 Exchange procedure 

Some of the samples are annealed in vacuum 

after polishing to remove the polising damage. 

The samples are placed in the Quarts tube and the 

system is pumped down to a pressure of 0.03torr. 

o o 

The samples are then annealed at 1100 Cl5 C for 30 
minutes. This also helps the sample and the 
Quarts tube to heat up and come into thermal 


68 



eauilibrium with the furn ace* For the 0-16 
anneal ? valve 2 is closed to cut off the vacuum 
pump and valve 3 is opened to introduce air into 
the system* So the old scale is grown in air at 1 
atmosphere pressure* The transfer of p-18 into 
the system for the tracer anneal is done without 
cooling down the system. At the end of 0-16 
anneal? valve 3 is closed and the system is 
evacuated with the vacuum pump by opening valve 2* 
Meanwhile the temperature of the zeolite tank is 
raised by using the heating coil* When the 
pressure of the system falls to 0.03 torr? valve 
2 is closed and the tracer is introduced into the 
system through valve 1. The 0-18 is kept at a 
slightly higher than atmospheric pressure to 
ensure no leakage of atmospheric air into the 
system which can cause dilution of the tracer* At 
the end of the tracer anneal? the system is 
evacuated again by cooling the zeolite tank with 
liauid nitrogen* On cooling it acts as a pump and 
sucks out all the tracer from the system. 


4.4 SIMS 



The GATAN SIMS at Pennsylvania State 
* 

University was used to obtain the data for the 
tracer concentration profile* The primary beam 
used was Ar* under a beam acceleration voltasfe of 
7KV. The complete instrument settings are listed 


below » 

Beam Acclerator 
Ion Source Filament 
Anode 

Beam Deflector Ini 

Lens Control Condensor 
Velocity filter 

Energy filter 


# Operated by Dr.C.Hauser 
Prof .C.Pantano. 


7 KV 1.53 mA 

17.2 A 
123 V 3.65 mA 

ial - 43 V Final - 204 V 
6.16 KV 

Electric Field- 0.81 V 
Magnetic field- 0.59 A 
Objective - 4.18 KV 

Outer potential- 40 V 
Inner potential- 113 V 
Bias - 39 V 

under the guidance of 
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The list of species simultaneously monitored and 


their mass to charge ratio is given below? 
Element Mass/Charge ratio monitered 


A1 

13*9 

16 

0 

16.4 

18 

0 

18.4 

Cr 

52.9 

Fe 

56.9 

Zr 

91.7 


Due to the insulating properties of alumina* 
it was necessary to compensate for the positive 
charge of the primary beam by spraying electrons 
on the sample surface* These electrons can cause 
desorpton of light elements from the near surface 
region* This gives rise to a background intensity 
of the monitored ions and is termed as Electron 

Stimulated Desorption(ESD) . The ESD level is 

/ 

evaluted by shutting off the primary beam and 
measuring the secondary ions due to the electron 
beam only. This background level is subtracted 
out from the actual reading. 
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A sample holder for SIMS can hold upto 6 

samples and has in addition 2 standards. The 

first standard is Ta 0 on Ta which is used to aim 

r> s 

the primary beam and to check if the rastering is 
correct. The second standard is ZnS which is used 
to aim the electron beam such that it is 
coincident with the primary beam. 

The samples are then sputtered one at a time 

and the output of the machine is fed directly into 

a computer. The output is in the form of the 

number of counts of a species versus the 

sputtering time. The conversion of counts to 

concentration of tracer C (ts) at a time t is 

18 s 

given by { 

h. R (t ) - 0-18esd 

c (t ) - l*—* ...(4.1) 

(h.„(t ) - 0-18ESD)+(h. At ) - 016ESD) 
lo s lo s 

where h (t ) and h (t > are the counts of 0-18 
18 s 16 8 

and 0-16 at a sputtering time t . 0-18ESD and 

s 

Q-16ESD are the ESD values of 0-18 and 0-16 
measured by shutting off the primary beam. The 
sputtering time is converted to sputtering depth 
by measuring the depth of the crater by a 
profilometer and dividing it by the total time 
under the primary beam giving the sputtering rate 
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S such thati 


♦ 

Y = S t ...(4.2) 

s s 

where Y is the sputtered depth eorrespondins to 
s 

the sputtering time t . 

s 
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CHAPTER 5 


RESULTS AND DISCUSSION 

5.1 Reasons for choosing the alloy and 
techniaue 

The theoretical modelling of the oxidation 

process of the FeCrAl-Zr alloy was done on the 

assumption that the oxide growth is planar? i.e.? 

the oxide metal interface is sharp. Smialek 

studied the microstructure of the oxide formed on 

o 

the same alloy. He found that at 1100 C? the oxide 

formed had closely packed polycrystalline grains 

and the oxide-metal interface was reasonably 

planar? at least for long oxidation times. At 
o 

1100 C there were very few pores observed in the 
oxide. Reddy C153 used the double oxidation 
techniaue on a NiCrAl alloy and concluded that the 
oxide was inward growing with the mechanism being 
oxygen grain boundary transport. There was no 
reason to beleive that the behaviour of FeCrAl-Zr 
alloy would be drastically different. So the 
FeCrAl-Zr alloy seemed like a good alloy to apply 
the double oxidation techniaue. Cawley C161 used 
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the double oxidation techniaue very successfully on 
Si and SiC. He used the SIMS to obtain the tracer 
concentration profile* He was able to 
experimentally obtain a sharp oxide interface as 
the oxygen signal in the SIMS fell down rapidly at 
the interface* He was also successful in 
determining the diffusion coefficients of oxygen 
from the observed concentration profile. His 
success was inspi rational in choosing the double 
oxidation techniaue and SIMS for these experiments* 

5*2 Lacey growth 

When the tracer concentration profile in the 

oxide on a FeCrAl-Zr sample was obtained* it 

differed from the profiles obtained by Reddy on his 

NiCrAl-Zr samples. To try to interpret the 

* 

difference the oxide was looked at under the SEM. 

The microstructure as shown in Figure 18 varied 

* 

drastically from that observed by Smialek on the 

same alloy. The sample had been polished by 3 

o 

micron diamond paster annealed in vacuum at 1100 C 

o 

for 30 minutes and oxidised for 40 hours at 1100 C. 

* Courtsey J. Smialek and J.Cawley at NASA Lewis and 
M.Occhionero at Case Western Reserve University. 
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FIGURE 18 - Top view of the oxide at the gas-oxide interface. 
Ohe sample has been polished to a 3 micron 
diamond paste finish, annealed in vacuum for 30 
minutes at 110093 followed by a 40 hour oxidation 
at 1100°C. 
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Instead of a compact polycrystalline sealer a lacey 


structure of about 2 micron wide ridses was 
observed at the oxide-gas interface. This shows 
that the oxide growth is not planar. 

To understand this difference in 

microstructure a systematic study was undertaken. 
The effect of three variables on the micro- 
structure was studied. the oxidation timer the 
surface finish of the sample and the time of vacuum 
anneal. The vacuum anneal is done before the 
oxidation to remove the polishing damage. 

The development of the oxide microstructure 

with time is shown in Figure 19. The alloy is 

polished to a 1 micron diamond paste finish and 

o 

oxidized at 1100 C for times of 3 hoursr 10 hoursr 
20 hours and 36 hours r 25 minutes. The lacey 
structure at the oxide-gas interface is evident at 
least for the shorter oxidation times. The ridges 
form a network of cells and the average cell size 
seems to decrease with time. The average cell size 
decreased from about 5 microns at 3 hours to about 
3 microns at 10 hours and goes down to about 1.5 
microns at 20 hours. At 36 hours and 25 minutes. 
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19 d) 36 hours,, 25 minutes oxidation at 1100°C. 

FIGURE 19 - Changes in microstructure of oxide at gas- 
oxide interface with increasing oxidation 
time. The samples have been polished to a 
1 micron diamond paste finish. 
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the lacea structure cannot be seen and the 
microstructure looks like that reported ba Smialek. 
So if onla the long time microstructure was 
examined it would seem to be a planar scale and 
would encourase the application of the model to 
derive Quantitative information. However* such a 
procedure would lead to errors as the oxide did riot 
start off looking that waa and at least the initial 
period of oxidation did not lead to a planar 
oxide growth. 

A detailed description of the oxidation 

mechanism leading to the lacea structure will not 

be attempted here. However it is important to note 

that the outward aluminum diffusion must occur to 

change the microstructure at the gas-oxide 

interface. So at least for oxidation up to 20 

hours* the oxide is not completela inward growing. 

Hindam and Smeltzer C491* observed lace on A1 0 

2 3 

scales grown on NiAl. Thea explained this 

occurence ba stating that the initial oxide formed 

was Y-Al 0 ba inward oxaSen diffusion. The scale 
2 3 

then went through a to tranformation accompanied ba 
a 14% reduction in volume. The a -alumina formed 
was for most part oriented. However* this volume 
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reduction Save rise to some faulted oxide regions. 
The scale was thus a combination of oriented 
alumina and misoriented oxide which provided the 
path for combined inward oxygen and outward 

aluminium diffusion# This gave rise to ridges at 
the gas - oxide and the oxide-metal interface 
correspondinS to misoriented oxide era ins. 

However? this does not explain why the cell size of 
the lacey structure should decrease with time as the 
t ransf ormat ion occurs only once. 

Hindam and Smeltzer observed that the 

structure of lace itself can be a function of the 
alloy grain orientation and can be vary drastically 
over two alloy grains in the same sample. Doychak 
E503 also observed a grain orientation dependency of 
the lacey structure during his study of NiAl 

oxidation. Work is presently being done by Smialek 
at NASA to understand the occurence of lace. The 
dependency of lace on the alloy grain orientation 
raises the Question whether the observed decrease 
in cell size with time is real or an artifact 
because the four samples were so chosen that the 
alloy grain orientation led to a decreasing cell 
size. Though this is a possibility? the odds 


83 



against it are high. 

Figure 19b) and Figure 20 show the effect of 
surface finish on the lacey structure* The sample 
in figure 20 has been polished by a 600 grit 
polishing paper while the sample in Figure 19b) has 
been polished by 1 micron diamond paste* Both the 

I 

samples have been oxidised for 10 hours at 1100 
C. No dramatic difference is observed in the lacey 
structure which suggests that surface finsh may not 
play a big role in lace formation. 

The effect of the vacuum anneal on the lacey 

structure can be seen by studying Figure 19b) and 

Figure 21* The sample in Figure 19b> has not been 

treated to the vacuum anneal while the sample in 

Figure 21 has been annealed in vacuum for 30 

o 

minutes at 1100 C. The sample with no vacuum 
anneal shows the lacey structure as a network of 
cells with the average cell size around 3 microns. 
However the structure of the sample with vacuum 
anneal looks drastically different. The lace takes 
up a spider like form with a nucleus and the lace 
spreading in radial directions from it* This shows 
up in certain areas of the oxide surface. Regions 
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FIGURE 20 - Top view of codde at the gas-oxide interface. 

The sample has been polised by 600 grit o polishing 
paper and annealed for 10 hours at 1100 C. 
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FIGURE 21 - Top view of oxide at gas-oxide interface. The 
sanple has been polished to a 1 micron diamond 
paste finish, annealed in vacuum for 30 minutes 
at 1100°C followed by a 10 hour oxidation at 
1100°C. 
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away from the lace covered areas have a needle like 

structure. A similar effect was observed at a 

higher temperature as shown in Figure 22. The 

sample in Figure 22a> was annealed for 30 minutes in 

o 

vacuum at 1200 C while the sample in Figure 22b) 

had not been given ana vacuum anneal. Both the 

samples were polished by a 1 micron diamond paste 

o 

and oxidised for 10 hours at 1200 C. The formation 
of lace on the sample with vacuum anneal was much 
more dramatic. These observations suggest that 
vacuum annealing encourages the formation of lace 
at least at short oxidation times. 

To study the effect of vacuum anneal on lacey 

structure for large oxidation times* two samples 

o 

were oxidised at 1100 C for 50 hours. Both samples 

o 

were annealed in vacuum for 30 minutes at 1100 C 
before oxidation. Figure 23 shows the 
microstructure of the oxide formed on the two 
samples. The sample in Figure 23a) has been 

polished with 1 micron diamond paste while the 
sample in Figure 23b) has been polished with 600 
grit polishing paper. Both samples do not show 
lace. However? Figure 18 shows a dramatically 
different microstructure. The sample in Figure 18 
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22 a) Vacuum anneal at 110 0°C far 30 minutes prior 
to oxidation. 
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22 b) No vacuum anneal. 


FIGURE 22 - Effect gf vacuum anneal on oxide microstructure 
at 1200°C. Both samples were polished to a 1 
micron diamond paste finish and oxidised for 10 
hours at 1200°C. 
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23 a) Microstructure of oxide at gas-oxide interface. The 
sample has been polished to a 1 micron diamond pagte 
finish, annealed in vacuum for 30 minutes at 1100°C 
and oxidised for 50 hours at 1100°C. 
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23 b) Microstructure of oxide at gas-oxide interface. Ihe 
sample has been golished by 600 grit polishing paper 
annealed at 110 0°C for 30 minutes and oxidised for 
50 hours at 1100°C. 


FIGURE 23 - Microstructure of oxide after long time oxidation. 
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has also been annealed in vacuum for 30 minutes at 
o 

1100 C after which it was oxidised for 40 hours at 
o 

1100 C. The specimen has been polished by 3 micron 
diamond paste. This difference in micro-structure 
maa be due to different alloy grain orientation. 

In conclusion! some of the facts about lace 
and its structure are stated without an attempt to 
explain the reason or the mechanism of structure 
variation. Occurence of lace seems to be sensitive 
to the vacuum anneal » at least at smaller times. 
The vacuum anneal encourages the formation of lace. 
The structure of the lace itself may vary from a 
cell like structure in Figure 14b) to a spider like 
structure ir» Fisure 21. The cell size and shape of 
the lace is sensitive to the alloy srain 
orientation. The cell size of the lace may 
decrease with time under isothermal oxidation. The 
lace seems to have a radial texture» i.e«> the 
ridses seem to form radially from a nucleus. 

The formation of lace and its structure does 
not seem to be consistent and can vary from sample 
to sample. Its erratic behaviour is a good reason 
to refrain from speculating on its origin. 
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5.3 Sputter etching 


An important point to consider is the aualita 
of data obtained ba SIMS. The Question to be 
addressed is that if the alloa was to undergo 
planar oxidation in a manner described ba the 
model t how accurate would be the tracer 
concentration profile obtained ba SIMS. Data 
obtained ba SIMS are in the form of ion counts 
versus sputtering time. Even though the counts are 
a strong function of the oxide matrix* the matrix 
dependenca is removed ba taking the ratio of 0-18 
counts to the total number of oxasen counts. 
Houever» the conversion of sputtering time to depth 
is valid onla if the sputtering rate is a constant 
or a known function of time. The manner in which 
data is collected ba SIMS is that the oxide is 
sputtered for a given length of time (1 sec in this 
case) and the counts reported is the average of the 
signal received over this period. 

To understand how SIMS collects the data* 
consider a planar oxide with a true concentration 
profile C(a) across it as shown in Figure 24. The 
data is collected ba sputtering for a timeAt and 
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FIGURE 24 - Sputtering of oxide scale by SIMS. 
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averaging the number of counts collected over that 
time period* Let A a » A« » ♦.««.£* be the 


12 n 

thickness of oxide sputtered away after successive 
sputterings of time At* The sputtered surface is 
assumed to be flat. Then the concentraton reported 
ba SIMS at time t» C(t>» is given baJ 


C(t)-C(^S±i) - i ( 


rr+1 


C dy 


...(5.1) 


If the sputtering rate is constant * then all the 
Aa's are the same. 


Figure 25a) shows a step function concentraton 
profile and Figure 25b> shows how the SIMS data 
would appear for a constant sputtering rate when 
the data in the form of counts versus sputtering 
time is converted to a plot of concentration 
versus distance. The profile obtained would be a 
sammetric distortion of the true profile. However 
it would be area preserving* i.e.* the area under 
the experimental concentration profile would be the 
same as that under the true concentration profile. 
Thus the mass balance ideas discussed in Chapter 2 
would be valid and’ applicable to the observed 
profile. However* if the sputtering rate is not 
uniform* then the conversion of counts versus time 
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to concent ration versus distance by usinsl Eauation 


4.2 would lead to an unsammet rical la distorted 
profile which would not be area preserving. In 
such a case the mass balance idea would not be 
applicable to the data* 

A further complication occurs if the sputtered 

surface of the oxide is not Planar but has instead 

a convoluted morphology as shown in Figure 26. Let 

C(a> be the true concentration profile of the 

tracer in the oxide. In this case for a sputtering 

time of At information is being collected from a 

range of depths y to y and the reported 

1 2 

concentration is the average of the concentration 

over this depth range weighted by the amount of 

material removed from a given depth. Let dV(y) be 

the volume of oxide removed at depth y by 

sputtering for a timeAt. The average concentrator! 

reported is given by? ^2 

_ / C(z) dV(z) 

C - 11 ...(5.2) 

I 2 

J dV(z) 
y l 

If the range of depths over which information is 
simultaneously received* is too large it may 
distort the concentration profile to an extent that 
obtaining Quantitative information from it is 
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volume of material dVfz) 



time At 


FIGURE 26 - Effect of non-uniform sputter surface. 
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impossible 


Another major problem faced when a sputter 
etchina techniaue is used on a polycrystalline oxide 
is the phenomenon of preferential sputtering. The 
grain boundary reslion may sputter at a faster rate 
than the grains. The sputtering rate of the grains 
themselves may be a function of their orientation. 
This means that the range of depths over which the 
signal is collected will vary across the sputter 
surface. Figure 27a) shows the oxide scale being 
made up of 2 regions A and B where B sputters twice 
as fast as A. A could represent the grain while B 
could represent the grain boundary or alternately 
they could represent grains of fast and slow 
sputter orientation. Figure 27 shows the sputter 
surface after a sputtering time Lt» The reported 
concentration is the average of the concentraton 
measured over the shaded area in Figure 27b). This 
pref erential sputter etching can also severely limit 
the resolution of the data. 

5.4 SIMS data 
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Figure 28 shows the total counts of oxygen 

plotted as a function of normalised scale depth for 

a sample polished by 1 micron diamond paste and 

oxidised successively in air and 0-18 tracer at 
o 

1100 C for 20 hours and 16 hours. 25 minutes 
respectively* It is important to note that the 
oxygen counts drop off very gradually and there is 
no discontinuity indicating a sharp oxide-metal 
interface* This could be due to a combination of 
preferential sputtering and non planar oxide 
growth. Figure 29 shows a SEM piciure of the same 
crater. The metal grains and the regions of oxide 
can be seen side by side as marked by M and 0 
respectively. The lacey structure is evident in 
the oxide. The apparent tracer concent ration 
profile of the same sample is plotted in Figure 30. 
If there was a sharp interface. the tracer 
concentration profile would have fallen off rapidly 
once the metal was reached. However, not only does 
the concentration profile at the interface not fall 
off steeply but it has a lesser slope than the near 
surface profile. This indicates that the problems 
of non planar oxide growth and preferential 
sputtering are too severe to apply our model to 
obtain Quantitative data. 
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0-16 oxidation time=20 hours 



sputtering time x 100 ( norma xi se d depth) 
total sputtering time 


FIGURE 28 - Plot of total oxygen counts versus normalised depth. 
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FIGURE 29 - Ibp view of the sputtered crater 





5.5 Conclusions 


The aim of this study was to develop the 
double oxidation tracer techniaue. The thrust of 
the work was to come up with a method of analysis 
of the tracer concentration profile such that 
valuable information could be obtained from it. 

The mass balance ideas developed can be used 
to deduce the oxidation mechanism from the tracer 
distribution in the oxide scale» if the oxide Srows 
in a planar fashion. For all planar growing 
scales, if a true tracer concentration profile can 
be obtained. the fraction of inward and outward 
growth can be deduced by mass balance. 

An effort was then made to develop an analysis 
which could be used to derive Quantitative 
diffusion data from the tracer concentration 
profile. The specific diffusion problem of tracer 
distribution in a polycrystalline oxide scale 
growing by inward oxygen grain boundary diffusion 
with exchange was salved using a model based on the 
Oishi-Ichimura combined lattice and boundary 
diffusion model. For uniform grains across the 


105 



scale* uniaue profles could be predicted using a 

combination of D* D and r. It was found that each 

b 

parameter had a different effect on the profile. 
The model could be generalized to accomodate any 
arain size distribution obtained in a real oxide 
scale. Careful consideration was also given to the 
factors that could distort the data. 

Having developed the analysis* it was tested 
on the FeCrAl-Zr alloy. Based on the extensive 
studies on NiCrAl alloys and long time oxidation 
studies on FeCrAl alloys by Smialek* the FeCrAl-Zr 
alloy was expected to oxidise in a planar fashion. 
However it was found that the oxidation behaviour 
of FeCrAl-Zr alloy was extremely unpredictable and 
the microstructure of the oxide varied from sample 
to sample and even within the sample. The 
formation and changes in the lacey structure at the 
gas-oxide interface suggested that outward cation 
transport also occured. So it was found that the 
FeCrAl-Zr alloy was not a good system to apply the 
model. The SIMS which had been very successfully 
used on amorphous materials led to preferential 
sputtering in the polycrystalline oxide scale and 
caused distortion of the tracer concentration 
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profile. Thus no attempt was made to derive 
diffusion coefficients from the data obtained. 

To exploit the analysis and Quantitative 
modellinsi done in this thesis? it is necessasary to 
apply the double oxidation techniaue to systems 
with typically planar growing scales. There is 

also a need to improve the profiling techniaues on 

1 

polycrystalline scales using proton activation or 
SIMS having uniform sputter etching. 
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CHAPTER 6 


SUGGESTIONS FOR FUTURE WORK 

1* The realistic oxide microstructure could be 
incorporated in the computer program by treating 
the problem of varying grain size and shape across 
the oxide scale# For oxide scales having columnar 
grains » the exchange process could be treated as 
diffusion into a cylinder# The only modification 
reauired would be to substitute Eauation 3.7 by the 
corresponding eauation of diffusion of tracer into 
a cylinder having a constant surface concent ration . 

2. The presence of precipitates of a second phase 
with different oxygen diffusivity could be 
incorporated in the program by considering the 
variation of lattice diffusivity across the oxide 
scale* 

3# The problem of combined inward oxygen tracer 
diffusion with exchange and outward aluminum 
diffusion could be addressed. 

4. The amount of 0-16 in the old oxide undergoing 
direct surface exchange with the gas during tracer 
oxidation can be calculated# This would lead to a 
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modification of the mass balance ideas. However 


the amount of 0-16 involved in direct exchange with 
the Sas is not very significant (Appendix B). 

5. A system with an inherently planar Srowins 
oxide scale should be chosen to apply this model. 
The analysis technicaues developed could be applied 
to the existing tracer concentration data on NiO by 
Atkinson et. al.C51U and NiCrAl-Zr by Reddy H15T. 
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APPENDIX A 


Listing of the Computer Program 


dimension cM0>»ci(40>»c3(40)»ti(40)»eSrain(40) 
open (unit^lFStatus^old' f name= ' fplf » dat ' > 

READ< 1 »*)concint 
read( 1 »*)doP»gc 
READ (If*) W3 

read( 1 1 * ) tstop » temp t tau » dl t gbwd* okp 

READ ( 1 i * ) U4 

read(l»#)gsize»dgb 

open (unit =2 t status* 'new ' r name= ' lpf 1 . ob J ' ) 
open ( uni t=3 » status= ' new '» name= ' grate . obJ 7 ) 
open ( unit=4 » status= 'new ' f name= ' gbc . ob J ' ) 
write(2»3) 

FORMAT < 5Xr 6 HTEMPf 1 10HD-18TIME t f 8H0~16TIME) 
write(2»4)tempftstopi't3u 
format(4XrF6«0»3:c»f6«0»4X»f6«0) 
write (2 » 8 ) 

FORMAT ( 5X * 1 4HGRA INSIZE < M ) f t 25HGB-DIFF 
C0EFF(M**2/SEC> f f 

21HLDIFF C0EFF(M**2/SEC> »5x»5hgbwd ) 
write(2f 9) ssize » dsb > dl t gbwd 

format <5x Fell .3»5x»ell.3» 13X»E11 .3? 10 k » el 1.3) 
write ( 2 t 10 ) 

f ornst < 3X 1 5H t ime 1 4;-: t Ihl * 5x t lh2 » 5x » 

lh4 1 5 x tlh6f 5:< t lhB, 5x f 2hl 0 f 4x f 

2hl2 » 4>: f 2h 14 f 4x f 2hl 6 f 4:: f 2h 1 8 f 4x f 2h20 ) 

READ (If#) U5 

read(l»*)tincFdtpr 

t=0, 

tprint=dtpr 
do 85 k = l > 20 
t i ( k ) = 0 . 0 
continue 
c( 1 )=0«5 
do 11 J=2 » 21 
c< J)=0.0 
continue 
do 700 J=1 » 20 
cgrain( J)=0.0 
continue 

tgrate=( <20.0/19.0>**2.0>*tau 

nscount=20 

xscsle=20 ♦ 0 




delty=( <okptftau>*>*0«f5} /.I.9 •> 0 
26 t=t+tinc 

3leha=ddb*dop/(sc*teniP*(<0kP*(T+t3u) >**0.5) > 
a=4*dl/(Ssize#dbwd> 

B=( (3.14/GSIZE)**2.)*DL*T 
var=0, 

DO 12 N4=l»20 
var=var+exp<-b# (N4#*2« ) ) 

VAR4=EXP(~B*(N4**2) ) 

IF ( 0AR4/VAR *LE ,0*01)G0 TO 100 
12 continue 

100 N4 = l 

beta=a*var 

>;=(dab/(delty**2.>>~(3lpha/(2.*de) ty) > 
y=( 1 , /tine)- (2*d3b/( del ty**2. ) )-beta 
z=(dsb/(delty#*2« ) ) + ( alpha/ (2 . *del ty ) ) 
do 16 J=2enscount 

16 ci ( J)=tinc*(;:*c< J+l >+y*c( J) 

% +z*c( J-l )+beta*cdrain( J> ) 

ci ( 1 ) =1 * 0*concir»t 
ci ( 40 ) =0 , 0 
do 35 k=l»nscount 
if ( ci (k ) » le* 0,0 )so to 35 
t i ( k ) =t i ( k ) + 1 i nc 
35 continue 

if (t-tprint)20f 17,17 

17 continue 

write (4»800)t>ci(2)*ci(8)»ci(16)r 

* ci(20)»c.t(21)rci<22>»ei(23>* 

* ci (24) » ci (25) >ci (26) »ci (27) »ci (28) 

800 format (lK»f7«0»12f5,2> 

do 80 J= i » 20 

b=( (3, 14/asize)#*2, )#dl#ti(J) 

call e;:pcal <b» vop) 

ca (J)=(l-( <6./C3.14**2 , ) )*vop) )* 

* (ci ( J)-carain( J) >+carain( J) 
ca( J) = ( ( (Ssize/2.0)#$3. )# 

* ca( J)+1.5*( (asize/2, >**2, >* 

* abwd*ci< J) )/( < (asize/2. )+(sbwd/2» ) >**3, ) 

80 continue 

write(2»18)t»ca(l>»C3<2>»ca(4)» 
t ca(6)»ca(8)»ca(10)>ca(12)sca(14)> 

* ca(16> »ca(18> »ca(20) 

18 format ( 1X»F7.0» 11F6*2) 
tprint=tprint+dtpr 

20 if ( t-tstop)515»21 >21 

21 do 90 J=l>r»scount 

b=( (3,14/asize)**2. )*dl#ti ( J) 
call expcal (b» vop) 
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* 

* 

* 

90 

* 

* 

30 

150 


600 


680 


* 

* 

610 

515 


530 


* 

* 

* 

580 

* 

* 

640 


ca< J)=(l-( (6./ <3. 14**2. ) )*vop> >* 

(ci(J)-c3rsin(J) )+c3rain(J) 

cs(J)=( ( <asize/2«0)**34 )# 

cs ( «j ) + 1 *5# ( (asize/2. )**2. >* 

abwd*ci ( J) )/( ( (asize/2. >+(abwd/2. > )**3. ) 

continue 

write(2f30)tfca(l)»c3(2)»ca(4)f 
ca(6>»ca(8)»ca(10)»ca(12)»C8<14)» 
ca ( 16 ) f ca ( 18 ) f CA ( 20 > 

FORMAT (1X»F7.0»11F6. 2) 

write(2»150)ci(20> 

format < ///f 7 « 2 ) 

areal=0«0 

do 600 n=lfl9 

areai= ( ca ( n ) +ca ( n+1 > ) /2 ♦ 0 

areal=are3l+3re3i 

continue 

s r e a 2 = 0 . 0 

nsc=nscount-l 

do 680 J=20»nsc 

sreaJ=l«0-( (ca(n)+ca(n+l) )/2.0) 

3re32=3rea2+3reaJ 

continue 

)<inc=(are3l-3res2)/(l40-ci (ns count) ) 

::scale=::5C3le+;;inc 

wri te ( 3 f 610 ) t f x scale t xinc t 

ca(21) »ca( 22) »ca< 23) »C3( 24) »C3< 25 ) f 

ca(26) rca(27) »ca(28) »ca(29) 

format ( lx f f 7. 0 f 2x f f 6 ♦ 2 f Ixf f 5 . 2 f lx f 9f 6« 2 ) 

stop 

continue 

tcounter=t+tau 

if (tcountei — ta rate ) 520 f 530 f 530 
continue 

do 580 J®1 mscount 

b=( (3« 14/2size)#*2« >*dl*ti(J) 

call execsl <bF vop) 

ca( J>=(1~( (6./(3.14**2« ) )*vop) )* 

(ci ( J)-c2rain( J) >+c2rain( J) 
ca( J)=( ( (3size/2«0>**3« >*ca( J) 

+1 ,5*( (Ssize/2« >**2. >* 

3bwd*ci( J) >/( ( (asize/2* )+(abwd/2. > )#*3. ) 
continue 

write (2f 640) tF cad ) >ca(2) ?ca(4) f 
ca(6) fC3<8) fc3 (10) Fca( 12) fcs( 14) 

Fca(16 ) f ca( 18 ) fC 3 (20) 
f ormat ( lx Ff7.0Fllf6.2) 
areal= : 0«0 
do 500 n=lFl9 
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areai=(c3(n)+c3(n+l ) )/2*0 
sres^aresl+areai 
500 continue 

sre32=0« 0 
Jnsc=nscount-l 
do 670 J=20»Jnsc 
aresJ=l .0-( (ca<n)+C3(n+l ) )/2.0) 
area2=area2+are3J 
670 continue 

;<inc=(3re3l-3res2)/( 1 ,0-ci (nscount ) ) 

xscale=xscale+xinc 

write <3*510)t rxscaler xinct 

* ca(21) »ca(22) » ca (23 ) »ca( 24 >»ca( 25 >f 

* c3(26)»ca(27)»ca(28)7ca(29) 

510 f ormat ( lx» f 7 . 0 > 2x f f 6 . 2 » lx ? f 5 . 2 r lx » 9f 6 . 2 ) 

tSrate= ( (tau**0.5)/19,0+(tSr3te**0.5> >**2,0 
cars in ( nscountt 1 >=ci (nscount > 
ci ( riscount+1 )=ci (nscount ) 
nscount =nscount+l 
520 ao to 24 

24 DO 25 K=1 » 40 
c(k)=ci (k) 

25 continue 
So to 26 
end 

subroutine expcsl ( bl > out ) 
out=0 « 0 
do 60 n=l » 20 
r.2=r.**2 

out =out+ ( exp ( -bl *n2 ) /n2 ) 
var2= ( out- ( e::p ( -bl #n2 ) /n2 ) > /out 
if ( var2.at »0.99)ao to 70 
60 continue 

70 r.=l 

end 
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APPENDIX B 


Justification of Assumptions Made 

While developing the mass balance ideas in 
Section 3.3> it is assumed that the amount of 0-16 
in the old scale involved in direct surface 
exchange with the tracer when the sample is 
oxidized in 0-18 gas is negligible. The diffusion 
of 0-16 atoms from the scale to the gas is due to 
the concentration Gradient only and can be 
evaluated by setting the chemical potential 
gradient eaual to zero. Area A in Figure 31 
represents the amount of tracer retained in the old 
oxide due to inward grain boundary diffusion of 
tracer with exchange. Area B is the amount of 0-16 
in the old scale which undergoes direct surface 
exchange with a same amount of 0-18 in the gas 
during the tracer oxidation. Area B is about 77. of 
area A. So the amount of tracer diffusing into the 
oxide due to direct surface exchange is not very 
significant . 
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FIGURE 31 - Comparison of direct surface exchange to total exchange. 




The problem of tracer exchange between the 


grain boundaries and grains was solved by treating 
the grains as spherical sinks with a constant 
surface concentration $ y (Section 3.5) ♦ This was 
done for mathematical convinience so that 4> y 
could be evaluated in £auation 3.8. However is 

not a constant but is a function of time. However 
the assumption of a constant grain surface tracer 
concentration is not very severe. As shown in 
Figure 32» the relatively high value of grain 
boundary diffusivity causes the grain surface 
tracer concentration to rapidly reach a near steady 
state value. There is only about an 8 % variation 
in the grain surface tracer concentration when the 
grains are being filled with tracer due to 
exchange. Thus the assumption of a constant 
does not lead to a large error. 
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FIGURE 32 


- The variation 


of grain surface and average grain 


concentration with time. 



16. Abstract 

Analysis techniques to obtain quantitative diffusion data from tracer concentration 
profiles were developed. Mass balance ideas were applied to determine the mecha- 
nism of oxide growth and to separate the fraction of inward and outward growth of 
oxide scales. The process of inward oxygen diffusion with exchange was theoret- 
ically modelled and the effect of lattice diffusivity, grain boundary diffusivity 
and grain size on the tracer concentration profile was studied. The development 
of the tracer concentration profile in a growing oxide scale was simulated. The 
double oxidation technique was applied to a FeCrAl-Zr alloy using 0-18 as a 
tracer. SIMS was used to obtain the tracer concentration profile. The formation 
of lacey oxide on the alloy was discussed. Careful consideration was given to the 
quality of data required to obtain quantitative information. 
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